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SUMMARY 



A method is described of rapidly estimating the posi- 
tion of the laminar separation point from the given pres- 
sure distribution along a body; the method is applicable 
to a fairly wide variety of cases. The laminar separation 
point is found by the von Karman-Millikan method for a 
series of velocity distributions along a flat plate, wh^ch 
consist of a region of uniform velocity followed by a re- 
gion of uniformly decreasing velocity. It is shown that 
such a velocity distribution can frequently replace the 
actual velocity distribution along a body insofar as the 
effects on laminar separation are concernod. 



An example of the application of the method is given 
by using it to calculate the position of the laminar sep- 
aration point on the N.A.C.A. 0012 airfoil section at zero 
lift. The agreement between the position of the separa- 
tion point calculated according to the present method and 
that found from more elaborate computations is "very good. 



INTRODUCTION 



Of the various available methods of calculating the 
characteristics of the laminar boundary layer in two- 
dimensional flow and, in particular, the position of the 
laminar separation point, the von Karman-Millikan method 
(reference 1) seems to be the most reliable. This method, 
when applied to the boundary- layer flow about an elliptic 
cylinder, showed good agreement with experiment when other 
methods failed (reference 2). The von K£ rman-Millikan 
method, however, has the disadvantage that the computa- 
tions are usually cumbersome and lengthy. Because of the 
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growing importance of the "boundary- layer prot>lom, soms 
method of rapidly finding the position of the laminar sep- 
aration point is needed. 

Although the calculations givon in reference 1 for a 
family of 11 doubl o- ro o f 11 velocity distribut-ions , consisting 
of a region in which the outside velocity increases lin- 
early with distance along the surface followed by a rogion 
in which the velocity decreases linearly, partly satisfy 
the need, some easy method having more goneral applicabil- 
ity is required. 

The purpose of the present paper is to develop a more 
general method of— rapidly estimating the position of the 
laminar separation point. The of-fectB of an adverse ve- 
locity gradient on laminar separation are studied. 

The method of estimating the position of the laminar 
separation point presented in this report is applicable to 
bodies having velocity distr ibut ions whose effects can be 
approximated by a region of uniform velocity followed by 
a region of uniformly decreasing velocity. It is to be 
not-ed that — the boundary- lay er velocity distribution at the 
point along the surface of a body where the outside veloc- 
ity is a maximum is very nearly the same as the Blasius 
distribution for a flat plate. Thus, the condit ion of the 
boundary layer at the point of maximum velocity can be re- 
produced by the flow over an equivalent length of a flat 
plate with 'uniform velocity equal to the maximum velocity. 
In a large number of cases, the velocity distribution back 
of the point — of maximum velocity can be well represented by 
a straight line up to the laminar separation point. The 
present method therefore has fairly general applicability. 

The laminar separation point is f~irst calculated by 
the von Karman-Mill ikan method for a series of velocity 
distributions over a flat plate; each distribution con- 
sists" 5f a region of uniform velocity followed by a- region 
of uniformly decreasing velocity. The velocity decrement 
sufficient to cause separation is then found a .a a function 
of the velocity gradient, which is expressed nondimens ion- 
ally. The decrement varies from 10.2 percent of the maxi- 
mum velocity to zero, depending on the value of— the gradi- 
ent . ! 

The relation of conditions on an airfoil section to 
those for which the calculations were carried out is dis- 
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cussed. As an example of the application of the method, 
the position of the laminar separation point on the N.A.C.A, 
0012 airfoil section at zero lift is calculated. 



CALCULATIONS AND RESULTS 



Throughout the present calculations, the reference 
length s 0 is taken as the distance from the leading edge 
of the plate to the point at which the adverse velocity- 
gradient is applied. The reference velocity U Q is taT-en 

as the velocity outside the boundary layer in the uniform- 
velocity region over the plate. 

The velocity gradient B 1 can then he expressed in the 
following nondimensional form: 

l = 5oW 

U 0 ds " - - 

where U is the velocity outside the boundary layer at 
any point along tho surface. 

s, tho distance along the surface. 

The following relations give the velocity distribu- 
tion outside the boundary layer over the plate: 

_tP- = 1 from = 0 to -— = 1 



and 



u o s o 
_f = i + , (ISLJo) . tor JL > i 



Figure 1 shows the form of the velocity distribution over 
tho plate for several values of 3P. 

The position of the laminar separation point was cal- 
culated by the method of reference 1 for a series of val- 
ues of P. In each case the velocity decrement Au/U 0 
sufficient to cause separation was found. The results of 
these calculations are given in figure 2 as a plot of 
AU/U 0 or U s /U 0 against P , where U s is the velocity 

outside the boundary layer at the separation point. 
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DISCUSSION... 



Figure 2 shows that small absolute values of the ve- 
locity gradient F correspond to large values of the 
velocity decrement AU/u o and that large values of F 
correspond to small values of Au/U 0 . When F is equal 
to zero, the outside velocity can be reduced by 10.2 per- 
cent of the maximum velocity U 0 before separation oc- 
curs, As the value of F increases,' the amount by which 
the velocity can be reduced after reaching its maximum 
valuo approaches zero. 

It is important to note that the magnitude of the di- 
mensional quantity dU/ds is not in it s elf- suf*f-ic ient to 
determine how much the vel-Qcity can be reduced from the 
maximum value before separation occurs. Figure 3 shows 
the effect of applying a_given velocity gradient du/ds 
at several distances from the leading edge. The' amount 
of the reduction from the maximum value before laminar 
separation occurs is much greater when the gradient Is 
applied close to the leading edge. .If the gradient _ la . 
applied immediately at the leading edge, the velocity dec- 
rement is independent of dU/ds. For this case F is 
always zero. 

When an attempt is made to apply the results of the 
foregoing calculations to an airfoil, the effects pro- 
duced by the actual velocity distribution over the body 
rrust be analyzed and compared with those produced -by the 
assumed type of velo city distribution. 

It is shown in reference 3 that, at any point along 
the surface of a body at which dp/ds is equal to zero 
(where p is the pressure), the curvature of the laminar 
boundary-layer profile at the surface is zero. This re- 
lation is applicable at the point of maximum velocity 
along a body as woll as to a flat plate with uniform ve- 
locity. Such considerations and comparisons with calcu- 
lations of actual boundary-layer profiles by the von 
Karman-Mi llikan method indicate that the shape ."of the 
boundary- layer profile at the point of maximum velocity 
outside the boundary layer is nearly the same as the 
Blasius flat-plate profile. Because any boundary- layer 

profile is specified completely by its shape an~d thick- 

ness, the effect of any region of rising velocity is only 
to affect the boundary-layer thickness at - the point of 
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maximum velocity. This effect can he reproduced hy a flow 
of uniform velocity over a suitable length o.f flat plate. 

The velocity distribution along the surface back of 
the point of maximum velocity can usually be approximated 
with sufficient accuracy by a straight line. The approxi- 
mation need be valid only as far as the point at which 
laminar separation occurs. 

Because a favorable pressure gradient makes for thin- 
ner boundary layers, the equivalent . length of flat plate 
is somewha.t less than the actual distance from the forward 
stagnation point of a body to the position of maximum ve- 
locity. The velocity over the plate is assumed equal to 
the peak velocity over the body. If boundary- lay er meas- 
urements are available at any point P near the point of 
maximum velocity, the length of plate equivalent to the 
region upstream of P is s x = §Rg/(5.53) a 

where 6 is the boundary- lay er thickness. 

Rg — US/i> , Reynolds Humher "based on the boundary- 
layer thickness. 

1?, the kinematic viscosity. 

The distance s 0 is then equal to s 1 plus the distance 
from P to the point s a at which the adverse velocity 
gradient is applied. (See fig. 4.) 

When no suitable boundary-layer measurements are 
available, the equivalent length of flat plate can be 
found from the following approximate relation: 

s/c 

c J \J 0 / VC/ 

0 

where U 0 is the peak velocity along the surface. 

c, the airfoil chord. 

The integration is to. be carried out over the region ex- 
tending from the forward stagnation point to the position 
of maximum velocity. In this case also, s Q = B 1 + s a . 
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The foregoing relation for the equivalent length of 
flat plate was derived from the assumption that the veloc- 
ity distribution in the boundary layer over the portion of 
the surface upstream of— the point of maximum velocity is 
similar to the Blasius flat-platre distribution and from 
the substitution of- this assumption in the von Earraan mo- 
mentum equation. " " " 

The integration indicated in equation (l) can, in gen- 
eral, be graphically performed. This graphical integration 
has been carried out for a series of velocity distributions 
similar to the distribution over the f-orward portion of the 
N.A.C.A.- 0.012 airfoil section at zero lift. for curves of 
this series, it was f"crund that" 



= 0.37 6 



•0.164 



3 m 



1m 



(2) 



where s m is the actual dis.tance along the surface 

from the forward stagnation point to the 
point of minimum pressure. 

p, the difference in pressure between tho points 



s n / 2 and s 



m* 



1m ' 



the local dynamic pressure at s n . 

Ap 



Equation (2) is valid for 0.003 < ~ < 0.2. 

Mil 

From the two straight lines representing the assumed 
velocity distribution, which replaces the actual velocity 
distribution over the body, the value of F is computr-ed: 



T = 



Sq dU 
U 0 ds 



o r 



1 = 



where p = (p - 



- s. 



dP 



, approximately 



2(1 - p m ) ds 

p co )/q is the pressure coefficient. 

p CT , the static pressure, i'n the undis- 
turbed stream. 
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p, the static pressure at any point 
on the airfoil. 



the dynamic pressure "based on the 
free-stream velocity. 

the pressure coefficient at the 
point of niniaun pressure. 

the slope of the line approximat- 
ing the pressure distribution 
hack of the point of minimum 
pressure. _ L 

The corresponding value of U s /Uo is fo -11 * 14 from figure 2. 
Separation is calculated to occur at the point where U 
falls to the value U 0 (U B /U 0 )'. 

As an example, the method has "been used to calculate 
the position of the laminar separat ion. po int on the H.A.C.A 
0012 airfoil section at zero lift. The velocity distribu- 
tion over the airfoil in terns of the free-stream velocity 
V was found by the method of reference 4. Figure 4 shows 
that, for unit chord, the position of the point of maximum 
velocity is 5j n = 0.125 and that Ap/q m = 0.037. Hence, 
from equation (2) , 



p m' 



dp 
ds ' 



s l _ - /_ v-0.164 



-±- = 0.376 (0.037) ' = 0.646 

and si = 0.0607. The position along the surface at which 

the adverse gradient is applied is s = 0.140. The dis- 
tance s 2 is then equal to 0.140 - 0.125 =0.015. The to- 
tal equivalent length of flat plate is b 0 = s x + s s = 
0.0957. The nondimensional velocity gradient la 

s Q dU 0.0957 , „ v 

J *%<u = " ri9o- (0 - 243) - - °- 0196 

Iron figure 2, U s /u 0 = 0.919, and the velocity at the 
separation point is 1,190 X 0.919 = 1.094. Separation is 
calculated to occur at s/c = 0.536, the position at which 
U/V falls to the value 1.094. — - — 
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The position of the laminar separation point was also 
computed according to the yon Karman-Millikan method. Two 
power series wore used to approximate the velocity distri- 
bution. These calculations indicated that the separation 
point was situated at s/c = 0.55. The agreement between 
the two sets of calculations is considered especially sat- 
isfactory In view of the simplicity of the present method. 

The method was also applied to .the case of the eLlip- 
tic cylinder report~BTL- i n reference 2. The previous calcu- 
lations showed that separation was to be expected at™ a 
distance s along the surface of 1.92. The present method 
gave s = 2.00 as the s eparat-ion point. The experimentally 
observed position was s = 1.99. 

CONCLUDING REMARKS 



A fairly general method of rapidly estimating the po- 
sition of the laminar separation point, based on the von. 
Kd.rma.n-Mi 11 ikan boundary- lay er theory, has been devised. 
Good agreement was obtain-ed between the results of this 
method and those of more elaborate computatrions when it 
was applied to caloulate; the laminar separation point on 
the N.A.C.A. 0012 airfoil section at zero lift. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , September 9, 1938. 
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